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Abstract 

We study the scalar glueball and vector meson spectral functions in 
a hot and dense medium by means of the soft-wall holographic model of 
QCD. Finite temperature and density effects are implemented through 
the AdS/RN metric. We analyse the behaviour of the hadron masses 
and widths in the (T, //) plane, and compare our results with the 
experimental ones and with other theoretical determinations. 

1 Introduction 

The properties of hadrons in a thermal and dense medium are currently 
widely investigated, both experimentally and theoretically. In this respect, 
the heavy-ion experiments at RHIC and LHC are producing a huge amount 
of data. However, from the theoretical point of view, it is still difficult to 
analyse these issues since one has to face nonperturbative calculations. 

The study of the spectral functions deserves particular attention since 
important information can be inferred from them, namely the behaviour of 
hadron masses and widths with temperature and density, or the highest tem- 
perature or density at which a state can exist. Moreover, from the vector 
meson spectral functions at low energy it is possible to compute transport 
coefficients using Kubo relations, together with other hydrodynamical quan- 
tities which can help in characterising the medium of quarks and gluons and 
in interpreting the experimental results ilij. 



1 



Hadronic spectral functions at finite quark chemical potential have been 
investigated by the Nambu-Jona Lasinio model [2l |3], QCD sum rules [H |5] 
and effective models [6H9]. On the other hand, this is a hard problem to 
deal with using lattice simulations, because of the sign problem occurring 
at nonzero chemical potential. Some analyses of vector mesons and scalar 
glueballs have been carried out by lattice simulations with two colours, where 
the sign problem is absent [T0Hl2]. A study with three colours in a mean 
field approximation at strong coupling has appeared in [13]. There are also 
experimental evidences of the modifications of p spectral function in medium 
[HHTH], as we discuss in Sect. |3| reviews on theoretical and experimental 
results can be found in [201 EI]- 

To address the nonperturbative problem of computing in-medium hadronic 
spectral functions, we use the AdS/QCD soft- wall model [22], which has 
been recently developed following the holographic bottom-up approaches to 
QCD inspired by the AdS/CFT (or gauge/gravity) correspondence [231425] . 
In these phenomenological approaches, gauge-invariant QCD operators are 
related to specific five- dimensional fields living in an anti-de Sitter (AdS) 
space, under the assumption that the strongly coupled gauge theory defined 
in a Minkowski space is dual to a semiclassical gravity theory in an AdS 
space. The duality consists in the equality of the partition functions of the 
two theories, under suitable boundary conditions, so that correlation func- 
tions of the strongly coupled theory can be obtained deriving a semiclassical 
functional [23]. This relation between a perturbative and a nonperturbative 
theory makes the correspondence interesting, in particular in view of the 
application to QCD. 

A key feature, that any holographic model aiming at describing QCD must 
incorporate, is the inclusion of a mass scale that breaks conformal invariance 
and is related to the scale of confinement. In the soft-wall model this is 
achieved by inserting in the action a dilaton background field containing a 
mass parameter ^22j. Other studies of hadronic spectral functions at finite 
chemical potential in so-called top-down holographic models can be found in 
[26| [27] . We shall focus on scalar glueballs and light vector mesons, but the 
method we adopt is quite general and can be applied to other states. 

2 Scalar glueball spectral function 

At T = and /i = scalar glueballs can be described in the soft-wall model 
by introducing a massless scalar field z), dual to the QCD operator 
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/3{as)G''^^G'''"' [281129], in the action: 

S = ^ld'xdz e-<^(^)y^^^^ (9mX) {d,,X) . (1) 

(j){z) is the dilaton term of the soft-wall model, and g is the determinant of 
the metric of the AdSs space 

ds' = ^ (dt^ - dx^ - dz^) , (2) 

with R the AdSs radius; in the following we set i? = 1. The holographic 
coordinate z runs in the range e < z < +oo, with e — )■ O"*". The dilaton 
profile is (j){z) = c^z"^, and the mass parameter c is fixed from the p meson 
mass predicted in the model: c = mp/2 = 0.388 GeV. The quadratic form of 
the profile has been chosen as the simplest one allowing to reproduce Regge 
behaviour for light hadron masses, and to get analytical results in several 
sectors, mostly at T = 0. The parameter k' is fixed from the matching with 
QCD of two-point function k' = n'^ / (IGa^Jf) (with /3i given in QCD by 
/3i = -llNjQ + Nf/3). 

In order to study the scalar glueball spectral function at finite tempera- 
ture and density, we follow the approach of [3T1433] , where the case of finite 
temperature and zero density have been considered; the results obtained and 
the tools developed can be generalised to the finite density case in a straight- 
forward way. 

When we switch temperature and quark chemical potential on, the action 
for the field X{x,z) is still ([T]), but with a different metric, which gathers 
information about temperature and density through a black-hole 



( dz^ \ 

ds^ = ^ \ f{z)dt^ - dx? - —— 0< z< Zh (3) 

f{z) = l-(^^ + q'zl^z' + qh' (4) 

where Zh is the position of the outer horizon of the black-hole, i.e. the lowest 
value which satisfies f{zh) = 0. This is known as AdS/Reissner-Nordstrom 
metric (AdS/RN). The temperature of the black-hole is 



47r 



df 



dz 



Z = Zh 



nzh 2 



being Q = qzf^^, with ^ Q ^ \/2. The parameter q represents the charge of 
the black-hole, and is related to the chemical potential p. Let us clarify briefly 



3 



this point. In the QCD generating functional, the quark chemical potential 
fi is the coefficient multiplying the quark number operator q'^q = qj^q. In 
the dual theory, such a coefficient is the source of the field associated to the 
quark number operator, namely the time component of a U{1) gauge field 
Am- The AdS/RN metric is the result of the interaction of this field with 
the AdS space. Instead of solving coupled differential equations involving the 
metric, the dilaton and the gauge field, as in [3ll |35] we use the ansatz: 

Aq{z) = fi — K q z'^ Ai{z)=0, i = 1,2,3, z, (6) 

and, from the condition AqIz^) = 0, we find 

H = KQ/zh. (7) 

K is a dimensionless parameter of the model, which must scale as a/A^ and 
can be determined studying different observables [35]; in the following we use 
K = 1. For convenience, we set c = 1, so all the following results are given in 
units of c. 

The bulk-to-boundary propagator K{p, z) of the glueball field in AD 
Fourier space is defined by X{p, z) = K{p, z)Xq{p), being X{p, z) and Xq{p) 
the bulk field and the source of the QCD operator, respectively. From the 
action ([T]) the equation of motion for K{p, z) can be computed 

~ A-f{z) + 2c'z^f{z) ~^,^ ( pI p'\~ 
K{p,z) K{p,z)^y—- — )^K{p,z) = ^; (8) 

as in |3T], the case p = will be considered, so uj^ = Vq- The solution must 
satisfy two boundary conditions: KiiJ^ , 0) = 1 and 

K{uj\ z)^{\- z/zn)"^^^ (1 + 0{l - z/zn)) , (9) 

i.e., near the black- hole horizon we select the in-falling solution. This latter 
condition is needed in order to obtain the retarded Green's function [361438] , 
which is computed by differentiating twice the on-shell action ([T]) with respect 
to the source A'o(w^). The result is [31] 



K Z 



(10) 



2 = 



Then, the spectral function is the imaginary part of the retarded Green's 
function 

p{uj') = 53 (ng(a;2)) . (11) 
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We have evaluated it for different values of temperature and chemical po- 
tential. Some numerical results are shown in Fig. [Tj in the left panel the 
first peak of the spectral function (multiplied by k') is plotted versus cu^ at 
T = 0.06c for different values of the chemical potential, whereas in the right 
panel the analogous curves computed at /i = 0.109c for various values of tem- 
perature are plotted. As already observed at zero chemical potential in [31] , 




Figure 1: Left: first peak of the scalar glueball spectral function at T = 0.06 
and fi = 0.109 (plain, blue), = 0.151 (dashed, red), /i = 0.191 (dotted, 
purple), = 0.221 (dash-dotted, green). Right: first peak of the spectral 
functions of the scalar glueball at /i = 0.109 and T = 0.06 (plain, blue), 
T = 0.07 (dashed, red), T = 0.08 (dotted, purple), T = 0.09 (dash-dotted, 
green). The variable p is defined as /) = k' p. All quantities are in units of c. 

at any finite and fixed value of p the peaks of the spectral functions broaden 
and move towards smaller values of cu^ as the temperature of the medium 
is increased, up to the melting point, where no peak can be distinguished 
anymore. The same behaviour is also observed at fixed temperature, increas- 
ing the chemical potential. The broadening is a sign that, as temperature 
and quark chemical potential increase, the states become unstable, getting a 
finite width. 

This effect can also be appreciated by a quantum-mechanical argument. 
Let us introduce the Regge- Wheeler tortoise coordinate r^, [33], defined by 
the relation dr, = —f{z)dz, ^ (r^, = 0, — oo correspond to 2; = 0, Zh, 
respectively). Then, through the Bogoliubov transformation X = e^^'^ip, 
with B{z) = z"^ + 3 log 2;, the equation of motion for the glueball field X 
becomes a Schrodinger equation for ip: 

- dl^in) + V{r,) ^(r,) = oo^tPir,) (12) 

with potential 

vi^)^{'^-'^)fu?-^mm^ (13) 
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The potential V is plotted in Fig. |2]as a function of r*. One can notice that 
at low T and fi, the potential has a well which makes the lowest eigenstates 
almost stable, while for higher values the decay probability becomes larger 
until the potential well disappears, as in [33]. 

4(1 k I I- 

30 - 
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Figure 2: Potential in the Schrodinger equation of motion for the scalar 



glueball, Eq.(12). The values of temperature and chemical potential are 
(c = 1): T = 0.02, /i = 0.1 (plain blue curve); T = 0.015, /i = 0.2 (dashed 
red curve); T = 0.06, /x = 0.1 (dotted green curve); T = 0.06, /i = 0.2 
(dash-dotted purple curve). 

In order to extract the dependence of glueball masses on the temperature 
and chemical potential, we have fitted each peak with a deformed Breit- 
Wigner function |311 132] : 

, , amV ,^ ^, 

Pbw[x) = 7 ^ , (14) 

where m is the mass of the state and F its width. We define the dissociation 
temperature as the lowest temperature at which the condition 



h , m F 

< 1 , h 



p2 ' (x — m?Y + m^F^ 



(15) 



is fulfilled. This means that we consider a state as melted if the ratio between 
the maximum of the pure Breit-Wigner function {h) and the squared width 
of the peak (F^), the two parameters characterising the form of the peak in 
the (m^,p) plane, is lower than 1. This criterion gives the same dissociation 
temperature at /i = found in [31], namely T ~ 0.116c. 
At very low values of both T and /z, we compute masses as eigenvalues of 
the equation of motion, since in this case the black-hole is far away along the 
holographic coordinate axis (z) and it does not affect the wavefunctions in an 
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appreciable way. We perform a different Bogoliubov transformation of tlie 
glueball field, z) = e^^'^/^H{p, z), with IS{z) = + 3\ogz- log f{z), 
such that the equation of motion for the field H{p, z) reads: 



-d'^,H{co\z) + U{z)H{co^,z) 
U{z) 



H{u\z) 



5/2 



B" 



(16) 
(17) 



The wavefunctions can be found by solving ([T6j), requiring that the eigen- 
function vanishes at 2; = and for large z. In Fig. [3] some of them are plotted 
at fixed T = 0.05c, for different values of the chemical potential. When /x is 
larger than a certain value, the solutions cannot be considered as eigenfunc- 
tions anymore, since their oscillatory behaviour near the black-hole horizon 
becomes important. We remark that the two methods described above to 
extract glueball masses give the same results for those values of T and for 
which both can be applied. 



H 




Figure 3: Glueball wavef unction at T = 0.05 and low values of the chemical 
potential: /i = 0.01 (blue plain curve), /i = 0.05 (red dashed), jj, = 0.09 
(yellow dotted), /x = 0.1 (green dot-dashed, all quantities in units of c). 
The latter is unstable, therefore for ^ 0.1 we cannot solve the eigenvalue 
problem to determine the glueball masses. 



The variation of the squared mass and width of the lightest glueball is 
shown in Fig. |4j we observe a decreasing of and an increasing of F with 
temperature and chemical potential. This behaviour is similar to what has 
been obtained by an analysis of scalar glueball masses in a medium in SU{2)c 
lattice field theory [12], in which a mass decrease at a very low temperature 
and high density in the transition from the hadronic to a superfiuid phase 
has been predicted. 
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Figure 4: Variation of the squared mass (left panel) and width (right panel) 
of the lightest scalar glueball with respect to temperature T and chemical 
potential fi. All quantities are in units of c. 

3 Vector meson spectral function 

An important application of the above described method is in the light vector 
meson sector. In fact, the p meson is considered as a very useful probe to 
investigate modifications of the properties of the hot and dense phase, since 
it can be easily reconstructed and its decay products carry information about 
the medium in the laboratory |20j . 

Following the soft-wall model formulation [22], we introduce the gauge 
fields Al and Aji associated to the SUl ® SUr gauge group, dual to the left- 
and right-handed quark currents of QCD, the global AD symmetries being 
gauged in the bulk. The vector field, which is dual to the QCD operator 
g7^g, is V = {Al + Aii)/2, and the action for this field in the quadratic 
approximation is 

S = -.rr^ [ d'xy^e-f(^^ Tr [F^'^FyMN] , (18) 

where F^^ = 9*^ - V^' , with V = 1/"T", T" being the generators of 
the SUv symmetry group; we shall choose the gauge Vz = 0. By comparing 
the two-point correlation functions of vector currents, computed in QCD 
and in the soft-wall model, one can fix kyg^ = 127r^/iVc [22]. The equation 
of motion for one of the three spatial components Vi in the Fourier space is 
(setting p = 0) 

f p-<f>{^) \ p-'t>{z) 

dz f{z)dMz,uj^) + -—^uj''V,{z,uj^) = 0. (19) 
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We compute also in this case the retarded Green's function, 



GfApo 



6^S 

6VP{-u)6V.^{uj) 



Sij 



-<t>{z) 



■ (20) 



z=0 



V{z, u"^) is the bulk-to-boundary propagator of the vector field in the Fourier 
space, defined by Vi{z,oj'^) = V{z,co'^) Vl^{u'^) and satisfying the equation of 



motion (19), with boundary conditions V{0,co^) = 1 and, as in the glueball 
case, the in-falling condition near the black-hole horizon: 



V(z,u') ^ (1 - z/z,)-'^i^^ (1 + 0{1 - z/z,)) . 



(21) 



Using Eq. (11), we find the spectral function in Fig. 5 for some values of 



temperature and chemical potential; the first peak corresponds to the p me- 
son. In particular the quantity p/u"^ has been plotted since we expect the 
vector spectral funtion to increase quadratically for large values of u (see, e.g. 
[39]). Like scalar glueballs, also vector mesons become unstable and finally 
melt while temperature and quark chemical potential are increased. In fact, 
at high enough T and/or p the spectral function becomes flat, with no more 
resonances, reflecting the formation of a weakly interacting plasma of quarks 
and gluons. Looking at Fig. [5| one can notice that the dissociation of the 
first excited state occurs much earlier than the one of the ground state (the 
p meson). 

The same analysis carried out in the glueball sector can be applied to 
vectors: using Eq. (14) it is possible to calculate the position of the peaks 
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Figure 5: Left: first two peaks of the vector meson spectral function at 
T = 0.06 and p = 0.139 (plain, blue), p = 0.225 (dashed, red), p = 0.259 
(dotted, purple), p = 0.455 (dot-dashed, green). Right: first two peaks of the 
spectral functions of vector mesons at p = 0.139 and T = 0.06 (plain, blue), 
T = 0.09 (dashed, red), T = 0.104 (dotted, purple), T = 0.15 (dot-dashed, 
green). All quantities are in units of c. 
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(masses) and their width, which are shown in Fig. [6j Again, at low values of 
the two parameters fi and T, the widths vanish and masses are computed by 
solving the equation of motion for the vector fields in the Schrodinger form. 




Figure 6: Variation of the squared mass (left panel) and width (right panel) 
of p meson versus temperature T and chemical potential /i. All quantities 
are in units of c. 

Fig. [6] shows that the p mass decreases with increasing temperature and 
chemical potential. In particular, at T ~ we observe a 13% effect (25% for 
the squared mass), while at an intermediate temperature T = 0.1c the effect 
reduces to 8% (16% for the squared mass). 

In the chiral model studied in ^40j to describe the in-medium behaviour 
of mesons, the scaling of the p-meson mass has been related to parameters 
governing chiral symmetry breaking. In particular, along with the dropping 
of the chiral condensate at finite density and temperature, a ~ 20% decrease 
of the p mass has been found at nuclear matter densities. In the soft-wall 
model both vector meson masses and the chiral condensate [H] decrease 
at increasing temperature and density, qualitatively in agreement with the 
Brown-Rho picture, suggesting the existence of a relation between these two 
phenomena. In our holographic model, chiral symmetry breaking and vector 
meson masses are two distinct sectors and a direct relationship probing such 
a scaling cannot be explicitly found: nevertheless, a comparison between the 
finite density behaviour of the pion decay constant and the p meson mass 
can be carried out, but deserves a separate study. 

In [42j, Hatsuda and Lee use QCD sum rules to obtain the density depen- 
dence of the p-meson mass; they find a decreasing behaviour fitted by a linear 
function: m(p)/m(0) = 1 — ap/po, where p is the density (po = 0.16 fm~^ 
is the nuclear matter density), and their prediction for a is a = 0.16 ± 0.06. 
Also in [l3] it is suggested that in-medium hadronic spectral functions carry 
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important information on chiral structure of hot/dense matter. By study- 
ing unitarized chiral models, a softening of a and p mesons is found, which 
is driven by the decrease of in nuclear matter. In [H] a quark model 
of nuclear matter is put forward, in which quarks interact with the nuclear 
field. As a result, a downward shift of 120 MeV has been predicted for the 
p mass. A mass dropping of p mesons has also been found in [7] within a 
coupled-channel analysis, in |8] within an extended vector meson dominance 
model, in |15] within sum rules at finite nucleon density, and in lattice SU(2) 
calculations [TOl [11] . However the behaviour of hadrons in medium is still a 
debated issue, both experimentally and theoretically. While the broadening 
of the spectral functions is widely accepted, the variation of the masses is 
unclear. Some data agree with a small decrease [151 [HI HE] , and other show 
a very small effect, even compatible with zero [T71 [19]. In particular, in |15] 
an analysis of the experimental data through various theoretical models has 
been carried out. In [16], a ~ 9% decrease in mass for the p meson at normal 
nuclear density is found from a linear fit, as suggested by the Hatsuda-Lee 
model [12] (a = 0.092 ±0.002). In [17] a modification of the p spectral func- 
tion is reported which is not compatible with a mass shift, but only with 
a broadening of the corresponding peak. In the analysis performed by the 
CLAS Collaboration [19] both the mass shifts by the Brown-Rho [40] and 
the Hatsuda-Lee [12] models are ruled out, since no evidence of large modi- 
fications to the p mass has been found (a = 0.02 ± 0.02), and an upper limit 
for a has been set to a = 0.053 with a 95% confidence level. 

The prediction of this model for the shift parameter a is a = (m(0) — 
m(/io))/m(0) ~ 0.012 at T = 0.023c = 9 MeV, where po = 0.209c = 
81.1 MeV is the value of the quark chemical potential corresponding to the 
nuclear matter density po- However, in the present model the quark den- 
sity p = K,q is related to the chemical potential in a model dependent way 
(through k), and again we have used k = 1. Although at values of p corre- 
sponding to the dissociation of the p meson we find a ~ 13% decrease of its 
mass, for p corresponding to nuclear matter density the effect is only ~ 1%, 
in agreement with the experimental analyses that find a small or compatible 
with zero effect. 

There are also theoretical studies predicting an increase in mass [G] [16] . 
A similar result has been found in the hard-wall model, another bottom- 
up holographic model characterised by a different mechanism of conformal 
symmetry breaking [17]; in particular, an increase of meson masses with 
density at T = is found by solving the corresponding equations of motion 
[^ . This behaviour, which is in contrast with the one found in our analysis 
and with the experimental results described above, is similar to the one found 
in other top-down holographic approaches [27] . 
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Finally, we remark that it would be interesting to extend the calculations 
based on the holographic approach to heavy mesons like J/V", investigating 
the behaviour of these states in a hot and dense medium; some methods have 
been proposed about how to describe charmonium-like states in the soft-wall 
model p|ll9]. 

4 Conclusions 

Studying scalar glueballs and light vector mesons in a hot and dense medium, 
we have observed a broadening of the peaks in their spectral function, which 
is a signal that the states become unstable as temperature and chemical 
potential are increased, and a shift of the position of the peaks towards 
smaller values of the mass, signaling a decrease of the hadron masses. While 
the first effect has also been found in other theoretical studies and in some 
experimental data, the latter is still a debated issue. We have also found 
that at T = the melting of glueballs and vectors occurs at ~ 0.365/€C 
and ~ 0.499kc, respectively. If we set k = 1/2, as in [31], the numerical 
values of the chemical potential are ~ 71 MeV for the glueball and fi ~ 
97 MeV for the vector mesons. At increasing temperature, these critical 
values decrease, and at T ~ 0.114c (T ~ 44 MeV) and T ~ 0.162c (T ~ 
63 MeV) the spectral function of glueballs and vectors, respectively, is almost 
flat already at zero density. Applying the soft-wall model, inspired by the 
AdS/QCD correspondence, we can address the nonperturbative problem with 
little effort, numerically solving a linear differential equation, and trying to 
collect as much information as such a phenomenological model of holographic 
QCD can provide. The broadening of the peaks confirms that the model 
captures an important accepted feature of QCD. The decreasing of the masses 
is a prediction which, particularly in the case of vector mesons, can be tested 
in the experiments; in the case of the density dependence, such a prediction, 
obtained in the soft-wall model, is different from the one obtained in other 
holographic approaches. 
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